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Abstract 

The  U.S.  Army  Corps  of  Engineers  Engineer  Research  and  Development 
Center  (USACE-ERDC)  Environmental  Lab  (EL)  assisted  USACE, 

Portland  District  (CENWP)  in  updating  a  CE-QUAL-W2  (W2)  model  of 
Applegate  Lake  based  on  a  previous  version  of  W2.  The  model  was 
calibrated  using  data  from  calendar  year  (CY)  2001  and  validated  with 
data  from  calendar  years  2003  and  2010.  One  set  of  W2  parameters  was 
successfully  applied  to  all  calendar  year  types  (2001  is  a  dry  year;  2003  is  a 
normal  year;  and  2010  is  a  wet  year).  This  model  and  the  corresponding 
results  from  the  study  provided  CENWP  with  more  refined  estimates  of 
water  temperatures  so  that  more  defendable  water  temperature  targets 
can  be  discussed  with  the  State  of  Oregon.  This  is  extremely  important 
because  the  Rogue  and  Applegate  Temperature  Total  Maximum  Daily 
Loads  and  Rogue  Spring  Chinook  Conservation  Plan  require  the  Corps  to 
review  the  Rogue  Basin  Project  operations  to  determine  whether 
improvements  can  be  achieved  to  downstream  temperature  for  the  benefit 
of  endangered  fish.  This  is  the  second  of  three  USACE  projects  on  the 
Rogue  River;  this  work  is  identical  to  the  Lost  Creek  Lake  Model  work  for 
CENWP. 
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1  Introduction 

Background 

Applegate  Lake  is  part  of  the  Rogue  River  Basin  Project.  Due  to  the  recent 
implementation  of  a  Biological  Opinion  on  the  system,  the  Rogue 
Temperature  Total  Maximum  Daily  Loads  (TMDL)  and  Rogue  Spring 
Chinook  and  Fall  Chinook  Conservation  Plans  require  that  all  systems  in 
the  Basin  be  reviewed  to  determine  whether  improvements  to  downstream 
temperature  targets  can  be  achieved  (ODEQ  20o8)(ODFW  2007HUSACE 
&  ODEQ  2009)  (ODFW  2013). 

Applegate  Lake  is  located  twenty-three  miles  southwest  of  Medford, 
Oregon,  on  the  Applegate  River  in  the  Rogue  River  National  Forest. 
Applegate  Dam  consists  of  an  earth-filled  embankment,  a  spillway,  a 
multi-level  withdrawal  tower,  a  regulating  outlet  conduit,  an  outlet  bypass, 
and  a  stilling  basin.  The  embankment  is  about  1,300  ft  long  and  about 
242  ft  high.  The  primary  authorized  purposes  of  the  dam  are  flood  control, 
fisheries  enhancement,  and  irrigation.  At  maximum  pool,  Applegate  Lake 
is  4.6  miles  long  and  stores  approximately  82,200  acre-ft  of  water  (USACE 
1990).  Figure  1  is  a  Google  Earth  screenshot  of  the  project  study  area. 

The  selective  withdrawal  water  temperature  control  tower  has  five  intake 
ports  that  allow  water  to  enter  one  of  two  wet  wells.  The  intake  inverts  are 
located  at  elevations  of  1,962  ft,  1,950  ft,  1,930  ft,  1,895  ft,  and  1,838  ft.  A 
diagram  of  the  multi-level  intake  tower  can  be  found  in  Figure  2. 

The  upstream  flows  into  the  lake  come  from  three  creeks:  Elliott  Creek, 
Middle  Fork,  and  Carberry  Creek.  The  reservoir  empties  in  to  the 
Applegate  River  near  RM  46.3. 

Approach 

In  order  to  determine  whether  new  temperature  targets  could  be  achieved, 
temperature  models  of  key  projects  needed  to  be  updated  and/or  created. 
Applegate  Lake  was  modeled  in  the  late  1980s  and  early  1990s  using  CE- 
QUAL-W2  (W2),  but  the  model  needed  to  be  updated  due  to  code 
revisions  and  operational  changes. 
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Objective 

The  goal  of  this  project  is  to  develop  and  calibrate  a  current  W2  model  for 
Applegate  Lake  that  can  also  be  used  to  fully  evaluate  the  effects  of 
operational  changes  on  release  temperatures  at  Applegate  Dam  on  the 
Rogue  River. 


Figure  1.  Google  Earth  image  of  the  Applegate  Lake  project  study  area. 


ERDC/EL  TR-17-6 


3 


Figure  2.  Representation  of  the  WTC  intake  structures  (Larson  1998). 
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2  Model  Selection  and  Development 

CE-QUAL-W2  (W2)  is  the  code  selected  to  develop  the  Applegate  Lake 
Model  (APPLM).  W2  is  a  2D  longitudinal-vertical  hydrodynamics  and 
water  quality  model.  It  is  capable  of  modeling  basic  eutrophication 
processes  and  is  best-suited  for  long,  narrow  waterbodies  that  do  not 
exhibit  substantial  lateral  variation.  W2  has  been  applied  to  hundreds  of 
studies  on  various  types  of  waterbodies  (rivers,  reservoirs,  lakes,  and 
estuaries)  all  over  the  world.  For  a  list  of  the  model  applications,  see  the 
CE-QUAL-W2  website:  http://www.ce.pdx.edu/w2/. 

CE-QUAL-W2  Description 

The  numerical  modeling  code  known  as  CE-QUAL-W2,  version  3.7,  was 
configured  for  application  to  Applegate  Lake.  W2  uses  a  finite  difference 
solution  of  the  laterally  averaged  equations  of  fluid  motion  (Cole  and  Wells 
2013).  It  allows  for  application  to  very  complex  water  systems  because  it 
accommodates  multiple  branches  and  multiple  waterbody  types.  W2  allows 
the  user  to  set  up  variable  grid  spacing  (longitudinally  and  vertically),  time- 
variable  boundary  conditions,  multiple  inflows  and  outflows,  and  time- 
variable  concentrations  for  each  water  quality  constituent  being  modeled. 
W2  (V3.7)  contains  a  user-defined  port  selection  algorithm,  which  allows 
the  user  to  specify  a  varying  number  of  elevations  for  dam  structures. 
Although  this  feature  is  not  utilized  in  the  calibration,  future  scenarios  may 
benefit.  In  addition  to  water  temperature,  W2  is  also  capable  of  modeling 
water  surface  elevation,  flow,  and  twenty-eight  water  quality  constituents 
such  as  total  dissolved  solids  (TDS),  inorganic  suspended  solids  (ISS), 
ammonium  (NH4),  biochemical  oxygen  demand  (BOD),  nitrate  (NO3), 
phytoplankton,  dissolved  oxygen  (DO),  and  organic  matter  (OM).  This 
study  only  modeled  temperature;  consequently,  the  other  constituents  will 
not  be  discussed. 

Project  Approach 

CE-QUAL-W2  is  well-suited  for  application  to  Applegate  Lake  for  the 
same  reasons  it  is  well-suited  for  application  to  Lost  Creek  Lake 
(Threadgill  et  al.  2015). 
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Three  in-lake  monitoring  stations  were  used  for  evaluating  model 
performance  during  calibration.  Locations  with  temperature  data  are: 
APGTi  (APP20001),  APGT2  (APP20002),  and  APGT3  (APP20003).  Data 
at  the  dam  and  downstream  from  the  dam  were  also  used  for  calibration. 
The  locations  of  the  sites  are  shown  in  Figure  3. 


Figure  3.  In-lake  profile  monitoring  stations.  Site  locations  provided  by  Kinsey  Friesen  (NWP). 


Calibration  Strategy 

Several  factors  were  used  to  determine  which  calendar  years  (CY)  were 
used  to  calibrate  and  validate  the  model.  The  largest  limiting  factor  was 
the  availability  of  observed  data.  Since  more  data  were  available  for  2001 
for  the  in-lake  stations,  CY01  was  used  to  develop  a  calibrated  model.  Once 
an  acceptable  set  of  calibration  parameters  were  found,  the  same  set  of 
model  parameters  was  used  for  CY03  and  CY10.  Each  of  the  chosen  years 
represents  various  water  year  types:  2001  was  a  dry  year;  2003  was  an 
average  year;  and  2010  was  a  wet  year. 
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3  Data  Analysis  and  Model  Preparation 

This  section  will  review  what  data  were  available  and  how  they  were  used 
to  define  the  calibration  input  files.  W2  has  several  data  requirements  to 
meet  before  simulations  can  begin: 

1.  Bathymetry  of  the  waterbody(ies) 

2.  Flow  and  temperature  characteristics  for  boundaries,  major  tributaries, 
and  point  sources 

3.  Dam  operations  and  structure  locations 

4.  Stage  data 

5.  Meteorological  conditions:  air  temperature,  dew  point  temperature,  wind 
speed,  wind  direction,  cloud  cover,  and  short  wave  solar  radiation  (if 
available) 

Model  Geometry 

Bathymetry  Data 

The  bathymetry  file  for  the  Applegate  Lake  Model  (APPLM)  was  originally 
developed  by  Mike  Schneider  (USACE)  for  the  original  W2  model  of 
Applegate  Lake.  Due  to  lack  of  available  documentation,  it  is  unknown 
where  he  obtained  the  bathymetry  data  (sediment  range  analysis,  cross 
sections,  etc.).  The  current  model  utilized  the  original  bathymetry  file, 
refined  the  grid,  and  modified  angles  of  the  original  segments  based  on 
Google  Earth  imagery. 

Model  Grid  Development 

Applegate  Lake  was  split  into  three  branches,  with  Branch  1  extending 
from  the  Applegate  Dam  at  the  Applegate  River  upstream  to  the  junction 
of  Middle  Fork  and  Elliott  Creeks;  Branch  2  constituting  a  side  channel 
that  enters  the  mainstem  of  the  reservoir  about  1.3  miles  upstream  from 
the  dam.  This  side  channel  has  inflows  from  Squaw  Creek,  but  due  to  the 
lack  of  data,  inflows  are  input  as  o  cms  in  the  model;  Branch  3  is  also  a 
side  channel  that  enters  the  reservoir  about  a  quarter  of  a  mile  upstream 
from  the  dam.  Its  main  inflow  is  the  French  Gulch,  but  again,  due  to  the 
lack  of  data,  this  inflow  is  also  input  at  o  cms.  The  reservoir  was  modeled 
with  41  longitudinal  segments,  varying  in  length  from  135.0  to  800.0  m, 
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and  77  vertical  segments  of  uniform  0.914  m  (~3  ft)  height.  Figure  5  is  a 
Google  Earth  image  with  the  segments  laid  out  as  the  model  sees  them. 

Table  1  provides  a  description  of  the  branches  in  the  reservoir;  the 
segment  numbers  do  not  include  the  inactive  (or  “null”)  segments  that 
start  and  end  each  branch  (required  in  W2).  Figure  4  shows  an  image  of 
the  longitudinal  segments  used  in  the  model  along  with  the  branch 
configuration.  Figure  5  is  a  Google  Earth  image  with  the  segments  laid  out 
as  the  model  sees  them. 


Table  1.  Geometry  characteristics. 


Description 

Branch 

Segment  Start 

Segment  End 

#  Segments 

Slope 

Branch  1  -  Mainstem  (Middle  Fork  and 
Elliott  to  Dam) 

1 

2 

29 

28 

0.000 

Branch  2  -  Ungauged  leg  of  the  lake 
(Squaw  Creek) 

2 

32 

35 

4 

0.000 

Branch  3  -  Ungauged  leg  of  the  lake 
(French  Gulch) 

3 

38 

40 

3 

0.000 

Figure  4.  Longitudinal  segments  with  branch 
configuration  for  the  APPLM. 
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Figure  5.  Google  Earth  image  with  model  grid  overlay  (produced  by  W2  Tools)  for  the  APPLM. 


The  bathymetry  of  the  APPLM  that  has  been  developed  has  been  verified 
to  replicate  the  observed  storage-elevation  curve  (obtained  from  NWP). 
The  storage-elevation  curve  obtained  from  NWP  was  dated  01/31/2006 
and  is  titled  "5%  Encroachment  on  Rule  Curve,  in  terms  of  Maximum 
Conservation  Pool.”  Figure  6  shows  the  storage-elevation  curve 
represented  by  the  model  compared  to  the  observed  storage-elevation 
curve  (or  volume-elevation  curve).  This  provides  the  ERDC  with 
confidence  that  the  bathymetry  is  good  and  sufficient  for  the  APPLM.  A 
complete  copy  of  the  bathymetry  file  used  can  be  found  in  Appendix  A  and 
the  model  input  file  was  delivered  to  CENWP. 

Dam  Features  and  Withdrawal  Locations 

Table  2  presents  an  abbreviated  list  of  segment  numbers  in  the  APPLM 
bathymetry  along  with  a  brief  description  of  what  site  is  located  at  the 
segment.  For  example,  the  in-lake  monitoring  site,  APPi,  is  located  at 
segment  28  in  the  APPLM  bathymetry. 
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Figure  6.  Volume-elevation  curve  comparison  for  the  APPLM. 


Table  2.  Model  segments  of  important  locations. 


Segment 

Length  (m) 

Distance 
Upstream  from 
Dam  (m) 

Distance 
Upstream  from 
Dam  (miles) 

Identification/Location 

1 

0 

7811.900 

0.000 

Boundary  (Null  Segment) 

2 

623.3 

7811.900 

6.621 

Beginning  of  Branch  1 

15 

298.7 

3274.400 

3.801 

In-lake  Station:  APP3  (APP20003) 

22 

274.3 

1503.800 

2.701 

In-lake  Station:  APP2  (APP20002) 

28 

159.1 

304.500 

1.956 

In-lake  Station:  APP1  (APP20001) 

29 

145.4 

145.400 

1.857 

DAM 

30 

0 

0.000 

1.767 

End  of  Branch  1 

31 

0 

2104.400 

1.767 

Beginning  of  Branch  2 

36 

0 

0.000 

0.459 

End  of  Branch  2 

37 

0 

739.100 

0.459 

Beginning  of  Branch  3 

41 

0.000 

0.000 

0.000 

End  of  Branch  3 
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Flow  and  Elevations 

Model  Inflow  Boundaries 

Upstream  and  Downstream  Boundaries 

Mean  daily  inflows  for  Applegate  Dam  were  available  only  for  the  year 
2003  and  later.  The  observed  inflows  were  obtained  by  NWP  from  the 
project  site  and  were  contained  in  the  gate  settings  file  spreadsheet.  For 
2001,  however,  these  data  were  not  available,  so  the  ERDC  obtained  six- 
hourly  estimated  inflow  via  the  NWPQuery  site  for  the  Applegate  Dam  site 
(APP).  W2  requires  that  all  branches  require  input  files  for  flow  and 
temperature.  However,  since  the  second  and  third  branches  are  ungauged, 
a  dummy  file  of  zero  flows  was  used  as  input  for  the  model.  These 
branches  were  included  in  the  model  only  to  capture  the  geometry  of  the 
reservoir  and  to  maintain  the  volume-elevation  relationship.  Their 
inclusion  will  have  no  impact  on  the  model.  The  model  will  fill  solely  using 
the  upstream  inflow.  Due  to  the  inaccuracy  associated  with  flow 
estimation,  a  decision  was  made  to  account  for  any  water  balance  issues 
(ungauged  flows,  rainfalls,  etc.)  by  using  the  water  balance  utility 
(available  with  the  W2  download). 

At  the  downstream  boundary,  located  at  the  dam,  total  outflows  were 
available  for  all  calendar  years  from  NWP  via  the  gate  settings  spreadsheet 
from  the  project  site.  Flow  at  each  intake  gate  was  recorded  each  day  as 
open,  closed,  or  the  amount  the  gate  was  open.  Based  on  that  information, 
the  gate  size,  the  total  outflow,  and  intake  gate  flow  was  calculated.  The 
elevation  data  available  at  the  dam  were  used  solely  for  model-to-data 
comparison.  Table  3  shows  the  data  sources  for  flow  and  elevation  for  the 
upstream  and  downstream  (Applegate  Dam)  boundaries.  Figure  7-Figure  9 
are  plots  of  all  flow  data  used  as  input  for  the  model  at  the  upstream  and 
downstream  boundary  for  all  three  calendar  years. 


Table  3.  Data  sources  for  flow  and  elevation  at  the  model  boundaries. 


River/Location  Name 

ID 

Source 

Variable 

Calendar  Year 

Upstream  Boundary 

NWP 

Flow,  Mean  Daily 

2001 

(Middle  Fork  and  Elliot) 

APP 

Gate  Settings  Spreadsheet  from 
Dam 

Flow,  Mean  Daily 

2003,  2010 

Downstream  Boundary 
(Applegate  Dam) 

APP 

Gate  Settings  Spreadsheet  from 
Dam 

Elevation,  Mean 

Daily 

2001,  2003, 

2010 

Downstream  Boundary 
(Applegate  Dam) 

APP 

Gate  Settings  Spreadsheet  from 
Dam 

Flow,  Mean  Daily 

2001,  2003, 

2010 
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Figure  7.  Flow  input  data  for  upstream  and  downstream  boundaries  for  CY01. 


Figure  8.  Flow  input  data  for  upstream  and  downstream  boundaries  for  CY03. 
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Figure  9.  Flow  input  data  for  upstream  and  downstream  boundaries  for  CY10. 


2010  -  Inflow  and  Total  Outflow  for  the  Applegate  Lake  Model 

Inflow  Outflow 


1  31  61  91  121  151  181  211  241  271  301  331  361 


Date,  MM/DD/YY 


Tributaries 

No  gauged  streams  discharge  to  Applegate  Lake.  For  this  reason,  no 
tributaries  were  defined  in  the  model.  Upon  initial  runs  for  all  calendar 
years  modeled,  the  model  still  seemed  to  underpredict  the  elevation  at  the 
dam.  For  this  reason,  a  distributed  tributary  was  added  to  the  system  to 
account  for  the  flow  imbalance. 

To  develop  a  distributed  tributary  input  file,  initial  model  output  and 
observed  elevations  must  be  input  into  the  Water  Balance  Utility  developed 
by  Portland  State  University  for  use  with  W2.  More  information  on 
developing  a  distributed  tributary  file  can  be  found  in  the  "Release  Notes” 
that  accompany  the  full  W2  download  along  with  the  Users’  Manual.  The 
water  balance  utility  will  calculate  negative  flows;  to  account  for  this,  the 
ERDC  averaged  out  the  negative  flows  with  the  positive  flows  from  the 
surrounding  time  period.  Figure  10  is  the  total  flow  that  was  added  to  the 
system  for  each  year  to  account  for  the  water  balance  problems. 
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Figure  10.  Distributed  tributary  Inflow  input  data. 


Distributed  Tributary  Inflows 

- 2001  - 2003  - 2010 


1 

A |  A  A 

£ 

353 


1-Jan  31-Jan  1-Mar  31-Mar  30-Apr  30-May  29-Jun  29-Jul  28-Aug  27-Sep  27-Oct  26-Nov  26-Dec 

Date,  DD-Mmm 


Model  Outflow  Boundaries 

The  amount  of  flow  withdrawn  through  each  intake  port  is  not  measured; 
however,  gate  settings  are  recorded.  Gate  setting  information  was  obtained 
from  NWP  as  an  Excel  spreadsheet  from  Dam  operators.  These  values 
were  then  used  to  develop  the  necessary  input  files  for  W2. 

Figure  11-Figure  13  are  plots  of  the  outflow  specified  at  each  intake 
structure  (intake  port).  The  ERDC  applied  conditions  to  the  total  outflow 
based  on  elevations  and  operations  procedures  as  detailed  in  the  Water 
Control  Manual  (U.S.  Army  Corps  of  Engineers  1991)  to  apportion  the 
total  outflow  to  each  intake  port. 
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Figure  11.  Outflow  input  data  at  specified  structure  for  CYOl. 


2001  Flow  into  Tower  as  Specified  in  Outflow  Input  File 


♦  STRl  •  STR2  ASTR3  ■  STR4  ♦  STR5  •  RO 


Date,  MM/DD/YY 


Figure  12.  Outflow  input  data  at  specified  structure  for  CY03. 


2003  Flow  into  Tower  as  Specified  in  Outflow  Input  File 


♦  STRl  •  STR2  ASTR3  ■  STR4  ♦  STR5  •  RO 


Date,  MM/DD/YY 
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Figure  13.  Outflow  input  data  at  specified  structure  for  CY10. 


2010  Flow  into  Tower  as  Specified  in  Outflow  Input  File 


♦  STRl  •  STR2  ASTR3  ■  STR4  ♦  STR5  *  RO 


Temperature 

Model  Boundaries 

For  all  calendar  years,  temperature  at  the  upstream  boundary  was  defined 
with  calculated  daily  inflows.  Temperature  was  not  measured  at  any  of  the 
upstream  sites  for  the  years  modeled  in  this  study.  The  ERDC  obtained 
period-of-record  daily  temperatures  at  the  Middle  Fork  Applegate  River 
(USGS 14361590),  Elliott  Creek  (USGS  14361600),  and  Carberry  Creek 
(USGS  14361700).  Data  were  available  for  10/01/1979-09/30/1987.  These 
mean  daily  temperatures  were  then  averaged  and  plotted  against  mean 
daily  air  temperature  from  Medford,  OR,  in  order  to  define  a  correlation  to 
use  to  estimate  inflow  temperature  for  the  modeled  years.  Figure  14  shows 
the  correlation  between  air  temperature  and  mean  water  temperature. 
Notice  the  R-squared  value  for  this  correlation  is  0.87,  which  suggests  that 
the  trendline  equation  (shown  in  the  chart)  represents  the  data  fairly  well. 
In  order  to  determine  just  how  well  the  equation  would  estimate  the  inflow 
temperature,  the  ERDC  used  the  correlation  to  estimate  what  the  tempera¬ 
ture  would  have  been  for  the  air  temperatures  in  1986;  Figure  15  shows  this 
comparison.  The  blue  line  represents  the  calculated  temperatures  using  the 
correlation  equation.  Overall,  ERDC  felt  the  equation  provided  a  good 
approximation  for  inflow  temperatures  when  no  data  were  available. 
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Figure  14.  Temperature  correlation  used  in  calculating  the  inflow  temperature. 


Figure  15.  Temperature  comparison  using  observed  and  correlated  temperatures  for  1986. 
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Temperature  at  the  upstream  boundary  was  also  used  as  input  for  the 
second  and  third  branches.  However,  since  flows  for  those  branches  are 
input  as  zero,  the  temperature  will  have  no  impact  on  the  model. 

Temperature  data  at  the  dam  were  used  as  calibration  data  for  the  model. 
Figure  16  provides  a  time-series  plot  of  temperature  at  the  upstream 
boundary  as  defined  in  the  model  for  all  calendar  years. 

Figure  16.  Temperature  input  data  for  the  upstream  boundary  for  2001,  2003,  and  2010. 
Estimated  Inflow  Temperature  for  Applegate  Reservoir 

- 2001  - 2003  - 2010 

1  31  61  91  121  151  181  211  241  271  301  331  361 

24  H - - - - - - - - - - - - h 


Date 


Tributaries 

Since  tributaries  were  not  monitored,  there  are  none  included  in  the 
model.  However,  because  a  distributed  tributary  must  be  used  to  improve 
the  water  balance,  the  upstream  temperature  input  file  was  duplicated  and 
used  as  input  temperature  for  the  distributed  tributary. 

Meteorological  Data 

The  same  meteorological  data  were  used  for  APPLM  and  APPLPM  that  was 
used  in  the  Lost  Creek  Lake  Model.  The  plots  of  the  data  are  shown  again 
below  as  a  reference.  Data  were  obtained  from  Medford,  OR.  Please  refer  to 
pages  16-19  °f  the  Lost  Creek  Lake  Model  Report  (Threadgill  et  al.  2015). 
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CE-QUAL-W2  Control  File 

The  control  file  for  the  model  calibration  (CYoi)  can  be  found  in  Appendix 
B  along  with  a  table  detailing  any  differences  for  all  other  model 
simulations.  In  order  to  keep  this  section  concise,  only  parameters  related 
to  temperature  will  be  discussed. 

Calculations,  Transport  Scheme,  and  Heat  Exchange 

Since  evaporation  is  always  considered  in  the  W2  surface  heat  exchange 
calculations,  it  is  important  to  turn  EVC  on  if  needed.  According  to  the 
manual,  if  calculated  inflows  are  used  in  setting  up  a  model,  then  EVC  is 
generally  set  to  OFF;  however,  in  the  case  of  the  APPLM,  EVC  is  set  to  ON 
since  we  are  using  direct  observed  inflows  from  the  project.  This  is  true  for 
all  three  years  modeled,  despite  the  fact  that  the  inflow  for  2001  was  an 
estimated  flow. 

The  transport  solution  scheme  used  in  the  APPLM  is  the  ULTIMATE 
scheme.  This  scheme  is  a  higher  order  solution  scheme  that  reduces 
numerical  diffusion  and  eliminates  the  over-  and  undershoots  that  the 
QUICKEST  scheme  generates  near  regions  of  shear  concentration 
gradients  (Cole  &  Wells  2013). 

In  the  W2  control  file,  the  user  must  specify  heat  exchange  parameters. 
The  first  parameter  specified  is  the  approach  used  for  computing  surface 
heat  exchange,  SLHTC.  For  the  APPLM,  the  ERDC  chose  to  use  SLHTC  = 
TERM  because  it  is  more  theoretically  sound  according  to  Cole  and  Wells 
(2013)  and  because  it  produced  better  model  results  than  SLHTC=ET. 
Since  the  meteorological  data  files  contain  shortwave  solar  radiation,  but 
the  model  produces  better  results  when  calculating  it  internally,  the  model 
setting  SROC  was  set  to  OFF,  which  specifies  that  W2  does  not  need  to 
read  an  extra  column  from  the  meteorological  input  file.  Although  the 
ERDC  was  provided  with  hourly  meteorological  data,  W2  was  still  allowed 
to  interpolate  the  input  data  to  correspond  to  the  model  time-step  by 
setting  the  parameter  METIC  to  ON.  The  wind  speed  measurement  height 
was  set  to  10  m  in  the  APPLM  as  indicated  by  the  14WS.  All  other  heat 
exchange  parameters  were  set  to  the  suggested  manual  values. 
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Extinction  Coefficients 

The  extinction  coefficient  card  contains  two  important  coefficients  for 
temperature  calibration.  When  water  quality  constituents,  other  than 
temperature,  are  not  being  modeled,  as  in  the  APPLM,  the  extinction  for 
pure  water,  EXH2O,  is  set  to  0.45  nr1  as  recommended  in  the  W2  manual. 
BETA,  the  fraction  of  incident  solar  radiation  absorbed  at  the  water 
surface,  is  also  set  to  the  value  of  0.45  in  the  APPLM  model.  The  W2 
manual  suggests  that  typical  values  for  BETA  are  approximately  0.2-0.7 
(Cole  &  Wells  2013). 

Selective  Withdrawal 

W2  is  capable  of  modeling  a  temperature  control  tower  with  selective 
withdrawal  features.  The  latest  version  also  has  the  added  capability  of 
dynamic  port  selection;  however,  since  this  was  not  used  for  the 
calibration  model,  it  will  not  be  discussed  here. 

The  Applegate  Lake  Water  Temperature  Control  tower  (WTC)  has  six 
intake  structures:  five  water  temperature  control  ports  and  one  regulating 
outlet  (RO)  (USACE  1990).  Figure  17  is  an  image  of  where  each  intake  port 
is  identified  in  the  model  control  file. 


Figure  17.  Schematic  representation  of  the  water  temperature  control  port  elevations. 
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4  Model  Calibration  Results  -  CY01 

Final  calibration  results  are  presented  in  this  section.  In  all  of  the  time 
series  plots  shown,  a  black  solid  line  represents  model  output,  a  solid  red 
circle  or  solid  or  dashed  red  line  represents  measured  data.  Three  statistics 
are  also  presented  in  the  charts:  mean  error  (ME),  absolute  mean  error 
(AME),  and  root  mean  square  error  (RMSE).  These  statistics  are  calculated 
as  shown  in  Equations  1-3.  The  model  was  output  every  day  as  a  daily 
average;  when  making  time  series  comparisons  to  the  observed  data,  a 
tolerance  of  0.5  days  was  used  for  the  model  output  so  that  model  output 
and  measured  data  were  compared  spatially  and  temporally  with  minimal 
averaging.  A  tolerance  of  7  days  was  used  for  the  model  output  when 
making  profile  plot  comparisons.  In  both  of  the  cases  with  the  tolerance  as 
selected,  the  statistical  comparison  is  a  one-to-one  comparison.  The  authors 
use  the  closest  date  and  the  closest  depth  for  comparing  values.  The 
tolerances  used  also  allowed  enough  spacing  to  avoid  observed  data 
averaging. 


y  (model  -  data) 

ME  =  — -  (1) 

n 

n 

Yj  abs(model  -  data) 

AME  =  — -  (2) 

n 


I  ( model  —  data)2 

RMSE  =  h  - - -  (3) 

Cumulative  distribution  plots  are  also  presented  in  this  section.  For  these 
plots,  the  solid  black  line  represents  model  output  and  the  dashed  red  line 
represents  observed  data.  These  plots  are  used  to  indicate  how  the  model 
is  behaving  overall  when  compared  to  the  observed  values.  For  example,  at 
high  temperatures,  the  model  over-/ underpredicts  temperature  by  XX 
deg-F,  where  XX  represents  the  AME  value.  Scatter  plots  are  also 
presented  to  give  a  statistical  representation  of  how  the  model  is  behaving. 
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A  general  rule  of  thumb  for  water  quality  calibration  is  that  the  absolute 
mean  error  should  be  within  10%  of  the  range  of  monitored  data  (Scott 
Wells)1,  temperature  AME  should  be  within  1  deg-C  (~i.8  deg-F),  and 
elevations  should  be  within  0.5  m  (1.64  ft).  Equation  4  is  the  equation  used 
to  calculate  the  target  values  for  AME.  These  target  values  were  calculated 
for  each  calendar  year  and  will  be  presented  in  tabular  form  in  the  following 
sections.  Units  for  these  targets  are  consistent  with  the  minimum  and 
maximum  values  for  each  constituent.  For  example,  for  flow,  the  minimum, 
maximum,  the  AME,  and  10%  target  are  presented  in  cubic  feet  per  second. 

Target  =  0.10* ((maximum  observed  value)  -  (minimum  observed  value))  (4) 


Flow 


Since  the  model  upstream  boundary  condition  segment  often  changes  based 
on  the  reservoir  volume,  the  ERDC  cannot  produce  flow  plots  to  verify  that 
the  upstream  boundary  condition  for  flow  is  satisfied.  Model  output  along 
with  observed  data  for  CY01  at  the  dam  is  shown  in  Figure  18.  Note  that  this 
is  really  just  a  representation  that  the  data  are  being  read  in  correctly  from 
the  input  outflow  file.  The  AME  for  all  data  pairs  for  2001  at  the  dam  is 
6.32  cfs,  which  is  well  less  than  1.0%  of  the  measured  range  of  flows  the 
calendar  year.  Table  4  presents  several  basic  stats  for  flow.  Based  on 
Figure  18,  the  slope  of  the  trendline  fitted  through  the  data  pairs  is  1.00  and 
the  R-squared  value  is  0.99.  Overall,  the  model  only  overpredicts  outflow  at 
the  dam  by  0.08  cfs. 


Table  4.  Basic  statistics  for  flow  (cfs)  for  CY01  calibration. 


SITE 

Observed 

Minimum 

Observed 

Maximum 

AME 

ME 

Slope 

R-Squared 

Dam 

51.00 

1210.00 

6.32 

0.08 

1.00 

0.99 

1  Wells,  Scott.  2008.  Personal  communication  with  Tammy  Threadgill.  June  15.  CE-QUAL-W2  Workshop, 
Portland,  OR 
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Figure  18.  Withdrawal  flow  at  the  dam  for  CY01  calibration. 


Temperature 

The  best  hope  in  correctly  predicting  the  outflow  temperature  is  to  correctly 
predict  the  in-lake  temperature  profiles  at  various  locations  in  the  reservoir. 
If  the  temperature  profiles  are  not  satisfactory,  the  chance  of  correctly 
predicting  total  outflow  temperature  is  highly  unlikely.  Profile  plots  and 
statistical  plots  for  all  in-lake  monitoring  sites  are  presented  in  Figure  19- 
Figure  24.  (Figure  3  shows  the  location  of  each  of  these  sites.)  A  time  series 
plot  and  statistical  plots  are  presented  for  the  dam  in  Figure  25.  The  average 
AME  for  each  of  the  in-lake  sites  are  within  the  acceptable  target.  Table  5 
presents  the  calculated  AME  and  the  temperature  target  that  ERDC 
attempted  to  reach  for  the  in-lake  sites  and  for  the  outflow  temperature  at 
the  dam.  Based  on  Figure  22-Figure  24,  the  average  slope  of  the  trendlines 
is  0.97  and  the  R-squared  value  is  0.95  for  the  in-lake  sites.  At  the  dam,  the 
AME  is  0.69  deg-C  with  a  slope  of  1.03  and  an  R-squared  value  of  0.99  (see 
Figure  25). 
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Table  5.  Basic  statistics  for  temperature  (deg-C)  for  CY01  calibration. 


SITE 

Observed 

Minimum 

Observed 

Maximum 

Target 

AME 

AME 

ME 

Slope 

R-Squared 

APP20001  (CYAVG) 

7.83 

24.54 

1.00 

0.80 

0.52 

0.88 

0.93 

APP20002  (CYAVG) 

10.29 

25.44 

1.00 

0.57 

-0.23 

0.95 

0.96 

APP20003  (CYAVG) 

19.14 

22.13 

1.00 

0.29 

0.19 

1.08 

0.96 

Dam  (Outflow) 

4.20 

19.80 

1.00 

0.70 

-0.52 

1.03 

0.99 

Figure  19.  Temperature  profiles  at  APP20001  in  CY01  calibration. 
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Figure  20.  Temperature  profiles  at  APP20002  in  CY01  calibration. 
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Figure  21.  Temperature  profiles  at  APP20003  in  CY01  calibration. 
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Figure  22.  Flow  linear  and  cumulative  distribution  plots  at  APP20001  for  CY01  calibration. 


Figure  23.  Flow  linear  and  cumulative  distribution  plots  at  APP20002  for  CY01  calibration. 


Figure  24.  Flow  linear  and  cumulative  distribution  plots  at  APP20003  for  CY01  calibration. 
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Figure  25.  Withdrawal  temperature  at  the  dam  for  CY01  calibration. 


Water  Surface  Elevation 

Model  output  along  with  observed  data  for  water  surface  elevations 
(ELWS)  in  CYoi  at  the  dam  is  shown  in  Figure  26.  The  AME  for  all  data 
pairs  for  2001  at  the  dam  is  0.64  ft  (-0.20  m).  Table  6  presents  the  basic 
statistics  for  observed  water  surface  elevation  at  the  dam.  The  slope  of  the 
trendline  fitted  through  the  data  pairs  is  1.01  and  the  R-squared  value  is 
1.0.  Overall,  the  model  only  overpredicts  ELWS  at  the  dam  by  0.55  ft. 


Table  6.  Basic  statistics  for  water  surface  elevations  (ft)  for  CYOI  calibration. 


SITE 

Observed 

Minimum 

Observed 

Maximum 

Target 

AME 

AME 

ME 

Slope 

R-Squared 

Dam 

1860.35 

1905.23 

4.49 

0.64 

0.55 

1.01 

1.00 
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Figure  26.  Water  surface  elevations  at  the  dam  for  CY01  calibration. 
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5  Calibration  Discussion 

Model  calibration  results  and  all  model  assumptions  are  discussed  in  this 
section.  As  stated  previously,  not  only  does  this  report  detail  graphical 
comparison,  but  the  authors  also  present  several  statistical  comparisons: 
AME,  RMSE,  and  ME.  Both  the  flow  results  and  the  temperature  results 
will  be  discussed  below. 

Water  Surface  Elevation 

As  stated  previously,  due  to  the  water  balance  instabilities  in  the  model,  a 
distributed  tributary  was  added  to  the  calibration  run.  This  drastically 
improved  the  initial  results.  Figure  27  shows  the  impact  of  not  using 
distributed  tributary.  Notice  how  the  model  severely  underestimates  the 
water  surface  elevation  for  ten  months  out  of  the  year.  By  the  end  of  the 
year,  the  model  needs  almost  8  ft  of  unaccounted  for  water.  Once  the 
distributed  tributary  was  added,  and  before  any  other  parameters  were 
modified,  there  was  definitely  an  improvement  to  the  results  (see 
Figure  28). 

Temperature 

Initially,  before  the  water  balance  issues  were  corrected,  the  model  was 
drastically  miscalculating  the  temperature.  However,  once  the  distributed 
tributary  was  added,  the  model  was  still  overpredicting  the  temperature 
(CYoi-Runo2).  Upon  observing  the  in-lake  profile  plots,  the  surface 
temperature  was  too  warm.  The  ERDC  performed  two  more  successive 
simulations  with  the  following  changes: 

1.  Decreased  the  shading  coefficient  from  1.0  to  0.85.  Due  to  the  fact  that  the 
lake  is  located  in  a  valley,  the  ERDC  has  decided  to  reduce  the  amount  of 
short  wave  radiation  actually  reaching  the  surface.  This  improved  the 
results  of  the  thermocline  location  in  the  profile  plots  at  the  in-lake  sites. 
(CYoi-Runo7) 

2.  Changed  EXH20  from  0.55  to  0.45  in  order  to  decrease  the  amount  of  heat 
retained  at  the  surface  instead  of  letting  the  heat  descend  into  the  water 
column.  Next,  the  team  changed  BETA  from  0.55  to  0.45.  BETA  is  similar 
to  EXH20  in  that  it  also  helps  to  retain  more  heat  surface.  These  changes 
had  a  very  small  positive  impact  on  model  temperature  predictions.  These 
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values  are  the  recommended  values  set  in  the  W2  manual;  however,  given 
the  fact  that  in  the  Lost  Creek  Lake  Model,  these  values  needed  to  be 
increased,  the  ERDC  originally  decided  to  leave  the  higher  values.  (CYoi- 
Runo8) 

Figure  27.  Time  series  and  statistical  plots  of  ELWS  without  the  distributed  tributary. 
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Figure  28:  Time  series  and  statistical  plots  of  ELWS  with  the  distributed  tributary. 


Temperature  comparisons  at  the  in-lake  stations  and  the  dam  between 
each  of  the  runs  discussed  above  is  seen  in  Figure  29-Figure  32.  In  all  of 
the  plots  below,  the  red  dots  are  observed  data.  The  time  series 
comparison  is  more  indicative  of  the  gains  in  temperature  improvement 
with  the  above  modifications  than  are  the  profile  comparisons. 


Figure  29.  Profile  comparison  at  APP20003. 
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Figure  30.  Profile  comparison  at  APP20002. 


Figure  31.  Profile  comparison  at  APP20001. 
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Figure  32.  Time  series  comparison  at  the  dam  for  CY01. 
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6  Model  Verification  Results  -  CY03  and 
CY10 

Model  verification  results  are  presented  in  this  section.  CY03  and  CY10 
were  used  because  they  had  the  same  types  of  monitored  data  and  similar 
available  in-lake  profile  data.  All  of  the  plots  and  statistics  presented  in 
this  section  were  developed  in  an  identical  manner  to  those  in  the  previous 
section. 


Flow 


Model  output  along  with  observed  data  for  CY03  and  CY10  at  the  dam  is 
shown  in  Figure  33  and  Figure  34.  Again,  this  is  really  just  a  representation 
that  the  data  are  being  read  correctly  from  the  input  outflow  file.  The  AME 
for  all  data  pairs  for  CY03  and  CY10  at  the  dam  is  121.05  cfs  and  124.65  cfs, 
respectively,  which  is  well  less  than  5%  of  the  measured  range  of  flows  for 
the  calendar  year.  Table  7  presents  the  5%  target  AME  and  the  model  versus 
observed  data  statistics. 

Temperature 

The  data  available  for  the  verification  years  was  a  little  different  than  in 
CY01.  For  CY03,  only  one  sample  date  at  all  stations  was  available 
(August  22).  For  CY10,  no  true  in-lake  stations  were  monitored.  In  order 
to  still  provide  feedback  on  in-lake  temperatures,  ERDC  chose  to  use 
temperatures  from  selected  dates  available  from  the  temperature  string 
located  at  the  dam  (in  place  since  2006).  It  is  important  to  note  that  the 
temperature  string  data  was  only  available  through  May.  The  15th  day  of 
Jan-May  was  chosen  as  representative  for  each  month  in  CY10.  Having 
plotted  the  observed  data  from  the  temperature  string  for  CY10,  the  plots 
showed  obvious  problems  with  the  gauge  at  elevations  1851  ft  and  1791  ft 
(see  Figure  39).  These  values  were  removed  from  the  plot  and  the 
statistical  comparison  (Figure  40). 


Table  7. 1%  Target  for  flow  (cfs)  for  CY03  verification. 


SITE 

Observed 

Minimum 

Observed 

Maximum 

Target 

AME 

AME 

ME 

Slope 

R-Squared 

Dam  -  2003 

99.00 

2520.00 

121.05 

25.44 

-2.77 

0.95 

0.95 

Dam  -  2010 

137.00 

2630.00 

124.65 

14.26 

3.07 

1.00 

0.99 
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Figure  33.  Withdrawal  flow  at  the  dam  for  CY03  verification. 


Figure  34.  Withdrawal  flow  at  the  dam  for  CY10  verification. 
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Profile  plots  and  statistical  plots  for  all  in-lake  monitoring  sites  are 
presented  in  Figure  35-Figure  41.  Time  series  plots  and  statistical  plots  are 
presented  for  the  dam  outflow  in  Figure  42  (CY03)  and  Figure  43  (CY10). 
Table  8  presents  the  calculated  AME  and  the  temperature  target  that  ERDC 
attempted  to  reach  along  with  comparison  statistics  for  the  in-lake  sites  and 
for  the  outflow  temperature  at  the  dam.  The  average  AME  for  all  of  the  in¬ 
lake  sites  are  within  the  acceptable  target  of  1  deg-C.  And  the  statistical 
values  are  within  typically  accepted  ranges.  At  the  outlet,  our  best  compari¬ 
son  site,  the  temperature  AME  is  0.80  deg-C  and  0.88  deg-C  for  CY03  and 
CY10,  respectively  (see  Figure  42  and  Figure  43).  The  model  underpredicts 
temperature  by  an  average  of  approximately  0.36  deg-C  at  the  dam. 


Table  8.  Temperature  stats  (deg-C)  for  verification  years. 


SITE 

Observed 

Minimum 

Observed 

Maximum 

Target 

AME 

AME 

ME 

SLOPE 

R-squared 

APP20003  (CY03) 

7.34 

23.91 

1.00 

0.37 

-0.01 

1.05 

0.99 

APP20002  (CY03) 

7.18 

24.12 

1.00 

0.53 

0.20 

1.03 

0.99 

APP20001  (CY03) 

7.02 

23.81 

1.00 

0.63 

-0.02 

1.05 

0.99 

Dam  Temp.  String  (CY10  AVG) 

4.83 

12.92 

1.00 

0.63 

-0.01 

1.20 

0.88 

Dam  -  Outflow  Temp  (CY03) 

4.86 

15.60 

1.00 

0.80 

-0.35 

1.09 

0.96 

Dam  -  Outflow  Temp  (CY10) 

4.50 

14.50 

1.00 

0.88 

-0.39 

1.02 

0.90 

Figure  35.  Temperature  profiles  at  in-lake  stations  in  CY03  verification. 
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Figure  36.  Flow  linear  and  cumulative  distribution  plots  at  APP20003  for  CY03 

verification. 


Figure  37.  Flow  linear  and  cumulative  distribution  plots  at  APP20002  for  CY03 

verification. 


Figure  38.  Flow  linear  and  cumulative  distribution  plots  at  APP20001  for  CY03 

verification. 
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Figure  39.  Temperature  profiles  at  the  dam  temperature  string  in  CY10  verification 

-  with  bad  observation  values. 


Figure  40.  Temperature  profiles  at  the  dam  temperature  string  in  CY10  verification. 
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Figure  41.  Flow  linear  and  cumulative  distribution  plots  at  the  dam  temperature 

string  for  CY10  verification. 


Figure  42.  Withdrawal  temperature  at  the  dam  for  CY03  verification. 
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Water  Surface  Elevation 

Model  output  along  with  observed  data  for  ELWS  CY03  at  the  dam  is 
shown  in  Figure  44  and  in  Figure  45  for  CY10.  Table  9  presents  several 
stats  and  lists  the  target  AME  for  each  verification  year.  In  general,  the 
ELWS  AME  should  be  within  0.5  m  (1.64  ft). 


Table  9.  Basic  statistics  water  surface  elevations  (ft)  for  CY03  verification. 


SITE 

Observed 

Minimum 

Observed 

Maximum 

Target 

AME 

AME 

ME 

Slope 

R-Squared 

Dam  (CY03) 

1882.89 

1987.06 

1.64 

0.84 

0.59 

1.00 

1.00 

Dam  (CY10) 

1882.18 

1986.97 

1.64 

1.10 

-0.07 

0.98 

1.00 
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Figure  44.  Water  surface  elevations  at  the  dam  for  CY03  verification. 
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Figure  45.  Water  surface  elevations  at  the  dam  for  CY10  verification. 
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7  Verification  Discussion 

This  section  serves  to  discuss  the  results  and  the  impacts  that  changes  have 
made  on  the  model  runs.  Just  as  with  the  Lost  Creek  Lake  Model,  no 
changes  were  made  to  the  control  file  for  either  verification  year.  Just  as  for 
CYoi,  a  distributed  tributary  was  needed  for  both  calendar  years.  As  stated 
previously  (see  Section  3,  Tributaries),  a  distributed  tributary  is  utilized  in 
W2  when  there  is  an  inconsistent  trend  with  the  water  balance  and  when 
the  user  can  account  for  missing  or  too  much  flow  (i.e.,  ungauged  flows).  It 
can  be  used  to  add  or  remove  water  from  the  system.  In  the  case  of  the 
APPLM,  a  distributed  tributary  was  used  to  add  water  to  the  system.  Refer 
to  Figure  10  for  the  flows  added  to  the  system  for  CY03  and  CY10. 
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8  Predictive  Port  Selection  Model 
Application 

In  order  to  provide  NWP  with  the  best  model  to  use  for  operation 
modifications,  the  calibrated  model  was  used  as  a  base  run  to  set  up  a  fully 
predictive  model  in  which  the  model  will  guide  dam  operations  based  on 
desired  temperature  targets.  The  temperature  target  presented  is  the  bi¬ 
weekly  target  developed  by  Oregon  Department  of  Fish  and  Wildlife  for 
2014  operations.  Based  on  results  from  the  Lost  Creek  Lake  Model 
(Threadgill  et  al.  2015),  the  ERDC-EL  chose  to  only  focus  on  the  USGS 
version  of  the  W2  model.  An  inventory  of  all  files  used  for  each  model 
simulation  can  be  found  in  Appendix  B  (Table  B3). 

USGS  -  W2  Predictive  Port  Selection 

Detailed  information  on  the  development  and  modifications  to  the  original 
W2  code  can  be  found  in  “Improved  Algorithms  in  the  CE-QUAL-W2 
Water-Quality  Model  for  Blending  Dam  Releases  to  Meet  Downstream 
Water-Temperature  Targets”  (Rounds  and  Buccola  2015).  Specifics 
relating  to  setup  of  the  Applegate  Lake  Predictive  Model  (APPLPM)  will  be 
discussed  here.  The  USGS  code  uses  an  iterative  process  to  determine  the 
optimal  flows  that  will  produce  the  desired  target  temperatures.  Of  course, 
this  means  that  the  run  time  will  also  increase.  In  the  case  of  the  APPLPM, 
using  this  code  tripled  the  run  time  (from  about  3-5  minutes  to  10-12 
minutes). 

Just  as  with  the  LCLPM,  There  were  no  changes  to  the  main  control  file 
from  the  calibration  model  (aside  from  output  filename  changes).  All 
completed  changes  were  made  in  the  w2_selective.npt  file,  which  is 
required  when  the  SELECTC  card  in  the  control  file  is  turned  ON. 
Although  the  structure  of  the  w2_selective.npt  file  is  very  similar  to  the 
PSU  version,  there  are  several  new  options.  The  new  cards,  not  discussed 
in  detail  here,  are: 

1.  NOUTS:  For  two  different  periods  throughout  the  year,  when  total  flow  is 
less  than  100  cfs,  the  RO  is  nonoperational.  To  specify  this,  NOUTS  during 
these  periods  was  set  to  5.  The  specified  period(s)  vary  between  different 
years.  The  time  intervals  used  for  each  year  can  be  seen  in  Table  10. 
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Table  10.  Intervals  when  R0  is  nonoperational. 


YEAR  -  TYPE 

TSTR 

TEND 

2001  -  DRY 

100.1 

335 

2003  -  AVERAGE 

32.1 

115.0 

325.1 

340.0 

2010  -  WET 

51.1 

74.0 

329.1 

346.0 

2.  TSSHARE:  For  the  APPLPM,  this  was  set  to  OFF.  DEPTH:  For  the 
APPLPM,  DEPTH  was  set  to  o  since  Applegate  Dam  consists  of  fixed 
ports. 

3.  MINFRAC:  For  the  APPLPM,  this  was  set  to  -2.832  cms  for  the  RO  only. 
The  negative  sign  indicates  that  this  an  actual  flow  value  instead  of  a 
percentage.  This  specification  is  made  due  to  the  minimum  gate  opening 
requirement  for  the  RO.  The  minimum  opening  is  0.6  ft;  although  the 
exact  flow  will  vary  based  on  elevation,  at  a  minimum,  flow  will  be  100  cfs 
at  the  lowest  elevations. 

4.  PRIORITY:  During  various  times  of  the  year,  NWP  operates  to  use  more 
surface  water  sometimes  and  at  other  times,  the  cold  lower  waters  are 
used.  Consequently,  for  the  fall  and  winter  months,  the  priority  was  shifted 
to  the  RO.  Outside  of  that  the  period,  priority  was  set  to  so  that  Intakes  1, 

3,  and  5  had  the  highest  priority,  and  Intakes  2, 4,  and  6  had  the  lowest 
priority. 

5.  MINHEAD:  For  the  APPLPM,  this  was  set  to  2.0  m.  There  is  about  a  6  ft 
minimum  head  according  to  the  WCM. 

6.  MAXHEAD:  For  the  APPLPM,  this  was  set  to  0.0  for  the  top  5  intakes.  For 
the  bottom  intake,  the  RO,  the  MAXHEAD  was  set  to  34.0  m.  This 
requires  the  elevation  to  be  a  conservation  pool  level  before  the  RO  is 
operational. 

7.  MAXFLOW:  For  the  APPLPM,  this  was  set  to  0.0. 

In  the  APPLPM  w2_selective.npt  file,  the  user  will  find  that  three  split 
times  were  identified.  The  reason  these  dates  were  identified  is  due  to 
operational  constraints  with  seasonal  withdrawal  depths.  Specifying  it  this 
way  allowed  ERDC-EL  to  set  the  PRIORITY  based  on  which  ports  were 
desired.  In  addition  to  the  w2_selective.npt  file,  because  DYNSEL  =  ON 
for  each  split,  the  user  is  required  to  have  an  additional  file  that  identifies 
the  target  temperatures  specified  by  Julian  day. 
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As  stated  above,  only  the  RO  had  a  MAXHEAD  value  specified.  Due  to 
operational  constraints,  the  RO  must  be  closed  or  open  no  less  than  0.6  ft, 
so  it  cannot  be  operated  unless  the  total  flow  is  greater  than  100  cfs  (2.832 
cms).  Even  with  the  additional  USGS  cards,  there  is  no  card  to  directly 
specify  this  requirement.  What  specifying  this  MAXHEAD  condition  does 
for  the  model  is  that  if  the  reservoir  is  more  than  34.0  m  above  the  RO,  the 
RO  will  not  be  used;  in  other  words,  the  RO  will  typically  be  used  when  the 
reservoir  is  at  the  lower  elevations.  Based  on  the  time  constraints  previously 
set,  the  RO  has  the  lowest  priority  during  the  spring  and  summer  already; 
so  setting  the  MAXHEAD  stipulation  also  adds  an  additional  check  on  the 
water  depth.  For  the  most  part,  at  Applegate  Lake,  the  RO  is  not  used 
frequently.  In  addition  to  this  MAXHEAD  constraint,  the  time  intervals 
during  which  the  RO  is  allowed  to  be  operational  varies  with  each  year. 

The  user  should  note  that  in  all  of  the  following  plots,  the  red  lines 
represent  a  temperature  target  range.  The  ODFW  targets  are  used  for 
determining  the  target;  however,  what  is  represented  on  the  following 
plots  is  a  target  range,  which  is  the  ODFW  temperature  target  +/- 1  deg-C, 
which  is  a  standard  measuring  error  for  temperature. 

Figure  46  and  Figure  47  is  the  w2_selective.npt  file  and  the 
dynsplit_selectiveX.npt  file,  respectively,  used  for  all  of  the  APPLPM 
model  runs.  In  Figure  48-Figure  80,  the  black  line/dots  represent  the 
results  from  the  calibration  run,  and  the  green  line/dots  represent  the 
results  from  the  predictive  mode  run.  Figure  48-Figure  54  are  plots  from 
CY01  (dry  year);  they  compare  the  results  from  the  calibration  and  the 
results  from  the  USGS-W2  blending  algorithm.  It  is  important  to  note  that 
at  no  time  in  the  calibration  model  were  Intakes  1-4  active;  during  the 
predictive  model  mode,  Intakes  4-6  were  active,  however.  Figure  55- 
Figure  67  represent  the  same  plots  for  CY03  (normal  year),  and  Figure  68- 
Figure  80  represent  CY10  (wet  year).  Figure  81  shows  the  average 
percentage  of  model-predicted  temperatures  that  fall  within  the  desired 
target  range.  As  one  can  see,  the  USGS  Port  Prediction  algorithm  produces 
better  results  more  often  than  the  calibration  using  the  gate  settings 
observed  flow  data. 
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Figure  46.  W2_Selective.NPT  file  used  for  the  APPLPM  for  CY10. 


blank  since  there  are  only  6  ports.) 


Figure  47.  dynsplit_selectiveX.npt  file  used  for  the  APPLMPM. 
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Figure  49.  CY01  -  Intake  4  -  temperature  into  tower. 
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Figure  50.  CY01  -  Intake  5  -  temperature  into  tower. 
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Figure  51.  CY01  -  RO  /  Intake  6  -  temperature  into  tower. 


Figure  52.  CY01  -  Intake  4  -  flow  into  tower. 


Figure  53.  CY01  -  Intake  5  -  flow  into  tower. 
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Figure  54.  CY01  -  RO  /  Intake  6  -  flow  into  tower. 
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Figure  55.  CY03  -  APPLPM  temperature  comparison  with  target  temperatures. 


Figure  56.  CY03  -  Intake  1  -  temperature  into  tower. 
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Figure  57.  CY03  -  Intake  2  -  temperature  into  tower. 


Figure  58.  CY03  -  Intake  3  -  temperature  into  tower. 


Figure  59.  CY03  -  Intake  4  -  temperature  into  tower. 
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Figure  60.  CY03  -  Intake  5  -  temperature  into  tower. 


Figure  61.  CY03  -  RO  /  Intake  6  -  temperature  into  tower. 


Figure  62.  CY03  -  Intake  1  -  flow  into  tower. 
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Figure  63.  CY03  -  Intake  2  -  flow  into  tower. 


Figure  64.  CY03  -  Intake  3  -  flow  into  tower. 


Figure  65.  CY03  -  Intake  4  -  flow  into  tower. 
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Figure  66.  CY03  -  Intake  5  -  flow  into  tower. 


Figure  67.  CY03  -  RO  /  Intake  6  -  flow  into  tower. 
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Figure  68.  CY10  -  APPLPM  temperature  comparison  with  target  temperatures. 
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Figure  69.  CY10  -  Intake  1  -  temperature  into  tower. 


Figure  70.  CY10  -  Intake  2  -  temperature  into  tower. 


Figure  71.  CY10  -  Intake  3  -  temperature  into  tower. 
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Figure  72.  CY10  -  Intake  4  -  temperature  into  tower. 


Figure  73.  CYIO  -  Intake  5  -  temperature  into  tower. 


Figure  74.  CYIO  -  RO  /  Intake  6  -  temperature  into  tower. 


ERDC/EL  TR-17-6 


58 


Figure  75.  CY10  -  Intake  1  -  flow  into  tower. 


Figure  76.  CYIO  -  Intake  2  -  flow  into  tower. 
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Figure  78.  CY10  -  Intake  4  -  flow  into  tower. 


Figure  79.  CYIO  -  Intake  5  -  flow  into  tower. 


Figure  80.  CYIO  -  RO  /  Intake  6  -  flow  into  tower. 
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Figure  81.  Average  %  of  model  temperature  within  the  target  range. 


%  within  Target  Range 


3 


Jan 

Feb 

Mar 

Apr 

May 

Jun 

Jul 

Aug 

Sep 

Oct 

Nov 

Dec 

■  W2-Calibration 

20.43 

40.48 

50.54 

0.00 

3.23 

35.56 

59.14 

16.13 

13.33 

36.56 

3.33 

11.83 

■  USGS-PredictiveRun 

20.43 

40.48 

51.61 

5.56 

34.41 

87.78 

100.00 

78.49 

63.33 

37.63 

0.00 

20.43 

ERDC/EL  TR-17-6 


61 


9  Summary  and  Conclusions 

The  USACE-ERDC-EL  assisted  CENWP  in  updating  a  W2  model  of 
Applegate  Lake  based  on  inputs  from  an  existing  model  of  the  reservoir. 
The  model  was  calibrated  and  verified  using  data  from  calendar  year  (CY) 
2001  (dry),  2003  (normal),  and  2010  (wet).  Across  all  calendar  years,  the 
model  captured  the  quantitative  and  qualitative  trends  for  temperature 
and  flow.  Quantitatively,  the  model  predicted  temperatures  within  1.0  deg- 
C  for  most  of  the  calibration  sites  (in-lake  sites  and  at  the  dam),  which  is 
far  better  than  many  other  temperature  studies  (Arhonditsis  and  Brett 
2004).  Qualitatively,  trends  were  consistent  with  measured  data.  Model 
performance  statistics  were  closely  paired  temporally  and  spatially  with 
the  measured  data. 

In  addition  to  a  fully  updated  calibrated  model,  the  ERDC-EL  also 
developed  an  application  of  the  model  using  modified  W2  code  from  the 
USGS  that  allows  for  a  better  functioning  blending  algorithm  between 
multiple  ports.  The  same  version  of  the  W2  model  was  also  used  for  a 
similar  application  of  Lost  Creek  Lake  (Threadgill  et  al.  2015). 

This  model  and  the  corresponding  results  from  the  study  provide  NWP 
with  a  fully  capable  model  that  helps  users  determine  how  operational 
changes  will  impact  downstream  water  temperature.  This  is  extremely 
important  because  the  Rogue  and  Applegate  Temperature  Total  Maximum 
Daily  Loads  (TMDL),  Rogue  Spring  Chinook  Conservation  Plan,  and 
possibly  the  Rogue  Fall  Chinook  Conservation  Plan  require  the  Corps  to 
review  operations  to  determine  whether  improvements  can  be  achieved  to 
downstream  temperature  for  the  benefit  of  endangered  fish. 

Additional  work  to  consider  is  a  deeper  investigation  of  the  actual 
attainability  of  the  target  temperatures  and  the  impacts  of  these 
temperatures  on  fish  with  respect  to  egg  emergence  data.  Based  on  the 
current  study,  there  will  be  times  during  the  year  when  reaching  the 
desired  temperature  target  is  simply  not  attainable  given  the  dam 
operation  criteria.  This  model,  coupled  with  an  in  depth  fish  analysis, 
would  provide  NWP  with  invaluable  information  regarding  dam 
operations  and  the  impacts  to  fish. 
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Appendix  A:  Bathymetry  File 

This  section  contains  an  image  of  the  bathymetry  file  used  for  the  APPLM. 
The  only  difference  between  calendar  years  was  the  initial  water  surface 
elevation  used  in  creating  the  bathymetry  file.  W2  V3.7  now  has  the 
capability  to  use  a  csv  file  developed  in  Excel.  The  images  below  (Figure  Ai- 
Figure  A3)  are  pages  from  the  Excel  file  used  to  develop  the  csv  file. 

Table  Ai  is  the  initial  water  surface  (ELWS)  used  in  the  development  of  the 
bathymetry  files  for  each  of  the  model  simulations. 


Table  Al.  Initial  ELWS  used  in  bathymetry  files  for  all  simulations. 


Calendar  Year 

ELWS  (m) 

ELWS  (ft) 

Calibration-2001 

569.585 

1868.72 

Verification-2003 

576.187 

1890.38 

Verification-2010 

577.044 

1893.19 
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Figure  Al.  Page  1  from  bathymetry  development  Excel  file. 
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Figure  A2.  Page  2  from  bathymetry  development  Excel  file. 
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Figure  A3.  Page  3  from  bathymetry  development  Excel  file. 
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Appendix  B:  W2  Control  File  with  Detailed 
Modifications  for  the  APPLM  and  the  APPLPM 

This  appendix  serves  to  present  the  control  file  (w2_con.npt)  used  for  the 
calibration  of  the  model  (see  Figure  Bi-Figure  Bio)  along  with  a  table  of 
changes  for  every  model  run  simulated  (see  Tables  B1-B3).  All  other  model 
simulations  will  be  compared  to  the  Calibration  w2_selective.npt  file  from 
CY01.  Discussions  of  all  modifications  are  made  in  the  main  report  text. 
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Figure  Bl.  Page  1  from  CY01  w2_con.npt  file. 


W2  Model  Version  3.7 

TITLE  C  . TITLE . 

Version  3.7  Applegate  Dam  Model 

Model  Run  from  January  1-Dec  31,  2001 

CY01-Run08  -  CY01-Run07  improved  surface  temp,  but  worsened  the  thermocline, 
therefore  decreasing  model  performance 

-  *  Thermocline  region  was  too  cool. 

-  **  Decrease  light  extinction  coefficients  to  allow  more  light 

to  penetrate. 

Tammy  Threadgill  -  USACE-ERDC-EL 

GRID 

NWB 

NBR 

IMX 

KMX 

NPROC 

CLOSEC 

1 

3 

41 

77 

1 

OFF 

IN/OUTFL 

NTR 

NST 

NIW 

NWD 

NGT 

NSP 

NPI 

NPU 

0 

6 

0 

0 

0 

0 

0 

0 

CONST ITU 

NGC 

NS  S 

NAL 

NEP 

NBOD 

NMC 

NZP 

0 

0 

0 

0 

0 

0 

0 

MISCELL 

NDAY 

SELECTC 

HABTATC 

ENVIRPC 

AERATE C 

INITUWL 

1000 

OFF 

OFF 

OFF 

OFF 

OFF 

TIME  CON 

TMSTRT 

TMEND 

YEAR 

1 

365 

2001 

DLT  CON 

NDT 

DLTMIN 

DLT I NTR 

1 

0.10000 

OFF 

DLT  DATE 

DLTD 

DLTD 

DLTD 

DLTD 

DLTD 

DLTD 

DLTD 

DLTD 

DLTD 

1.00000 

DLT  MAX 

DLTMAX 

DLTMAX 

DLTMAX 

DLTMAX 

DLTMAX 

DLTMAX 

DLTMAX 

DLTMAX 

DLTMAX 

3600.00 

DLT  FRN 

DLTF 

DLTF 

DLTF 

DLTF 

DLTF 

DLTF 

DLTF 

DLTF 

DLTF 

0.90 

DLT  LIMI 

VISC 

CELC 

wb  1 

ON 

ON 

BRANCH  G 

US 

DS 

UHS 

DHS 

UQB 

DQB 

NLMIN 

SLOPE 

SLOPEC 

BR1 

2 

29 

0 

0 

0 

0 

1 

0.00000 

0.00000 

BR2 

32 

35 

0 

19 

0 

0 

1 

0.00000 

0.00000 

BR3 

38 

40 

0 

27 

0 

0 

1 

0.00000 

0.00000 

LOCATION 

LAT 

LONG 

EBOT 

BS 

BE 

JBDN 

WB  1 

42.0565 

123.115 

539.667 

1 

3 

1 

INIT  CND 

T2I 

ICEI 

WTYPEC 

GRIDC 

WB  1 

4.90 

0.00000 

FRESH 

RECT 

CALCULAT 

VBC 

EBC 

MBC 

PQC 

EVC 

PRC 

WB  1 

OFF 

OFF 

OFF 

ON 

ON 

OFF 

DEAD  SEA 

WINDC 

QINC 

QOUTC 

HEATC 

WB  1 

ON 

ON 

ON 

ON 

INTERPOL 

QINIC 

DTRIC 

HDIC 

BR1 

ON 

ON 

ON 

BR2 

ON 

ON 

ON 

BR3 

ON 

ON 

ON 

HEAT  EXCH  SLHTC 

SROC 

RHEVAP 

METIC 

FETCHC 

AFW 

BFW 

CFW 

WINDH 

WB  1 

TERM 

OFF 

OFF 

ON 

OFF 

9.20000 

0.46000 

2.00000 

10.0000 

ICE  COVE 

ICEC 

SLICEC 

ALBEDO 

HWICE 

BICE 

GICE 

ICEMIN 

ICET2 

WB  1 

OFF 

DETAIL 

0.25000 

10.0000 

0.60000 

0.07000 

0.05000 

3.00000 

TRANS POR 

SLTRC 

THETA 

IWB  1  ULTIMATE 

0.55 

HYD  COEF 

AX 

DX 

CBHE 

TSED 

FI 

TSEDF 

FRICC 

ZO 

WB  1 

1.00000 

1.00000 

0.30000 

11.984 

0.01000 

1.00000 

MANN 

0.00100 

EDDY  VISC 

AZC 

AZSLC 

AZMAX 

FBC 

E 

ARODI 

STRCKLR 

BOUND FR 

TKECAL 

WB  1 

W2 

IMP 

1.000 

N  STRUC 

NSTR 

BRl 

6 

BR2 

0 

BR3 

0 

STR  INT 

STRIC 

STRIC 

STRIC 

STRIC 

STRIC 

STRIC 

STRIC 

STRIC 

STRIC 

BR  1 

ON 

ON 

ON 

ON 

ON 

ON 

BR  2 

BR  3 

STR  TOP 

KTSTR 

KTSTR 

KTSTR 

KTSTR 

KTSTR 

KTSTR 

KTSTR 

KTSTR 

KTSTR 

ERDC/EL  TR-17-6 


69 


Figure  B2.  Page  2  from  CYOl  w2_con.npt  file. 
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Figure  B3.  Page  3  from  CY01  w2_con.npt  file. 
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Figure  B4.  Page  4  from  CY01  w2_con.npt  file. 
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Figure  B5.  Page  5  from  CY01  w2_con.npt  file. 
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Figure  B6.  Page  6  from  CY01  w2_con.npt  file. 
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Figure  B7.  Page  7  from  CY01  w2_con.npt  file. 
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Figure  B8.  Page  8  from  CY01  w2_con.npt  file. 
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Figure  B9.  Page  9  from  CY01  w2_con.npt  file. 
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Figure  BIO.  Page  10  from  CY01  w2_con.npt  file. 
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Table  Bl.  Changes  to  Calibration  w2_con.npt  File  for  Other  Runs. 


RUN 

YEAR 

TEMPI 

TSED 

Calibration-2001 

2001 

4.90 

11.984 

Verification-2003 

2003 

6.20 

12.513 

Verification-2010 

2010 

4.70 

11.743 
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Table  B2.  Inventory  of  files  needed  to  run  the  LCLM. 


Run  Name 

CY01_Run08 

CY03-Run02 

CY10-Run02 

File  Type 

Calibration  - 
2001 

Date  Stamp 

Verification  - 
2003 

Date  Stamp 

Verification  - 
2010 

Date  Stamp 

W2_C0N.NPT 

“ 

10/26/15 
2:10  PM 

“ 

01/09/15 
10:57  AM 

“ 

10/26/15 
3:57  pm 

GRAPH. NPT 

- 

10/22/12 
5:01  PM 

- 

10/22/12 
5:01  PM 

- 

10/22/12 
5:01  PM 

WSC  File 

APP-WSC.NPT 

12/23/14 
2:30  PM 

APP-WSC.NPT 

12/23/14 
2:30  PM 

APP-WSC.NPT 

12/23/14 
2:30  PM 

SHD  File 

APP-SHD-2.NPT 

01/05/15 
10:49  AM 

APP-SHD-2.NPT 

01/05/15 
10:49  AM 

APP-SHD-2.NPT 

01/05/15 
10:49  AM 

BTH  File 

APP-BATH-FINAL- 
2001. NPT 

12/04/15 
2:21  PM 

APP-BATH- 
FINAL- 
2003. NPT 

01/14/15 
12:41  PM 

APP-BATH- 
FINAL- 
2010. NPT 

01/08/15 
1:06  PM 

MET  File 

APP-MET- 
2001. NPT 

01/27/14 
10:53  AM 

APP-MET- 
2003. NPT 

02/03/14 
2:06  PM 

APP-MET- 
2010. NPT 

01/14/15 
2:54  PM 

QIN  File 

APP-QIN- 
2001. NPT 

12/04/14 
1:44  PM 

APP-QIN- 
fromGates- 
2003. NPT 

03/04/14 
2:29  PM 

APP-QIN- 
fromGates- 
2010. NPT 

01/14/15 
2:49  PM 

APP-BR2- 

QIN.NPT 

12/17/12 
4:18  PM 

APP-BR2- 

QIN.NPT 

12/17/12 
4:18  PM 

APP-BR2- 

QIN.NPT 

12/17/12 
4:18  PM 

APP-BR3- 

QIN.NPT 

12/17/12 
4:18  PM 

APP-BR3- 

QIN.NPT 

12/17/12 
4:18  PM 

APP-BR3- 

QIN.NPT 

12/17/12 
4:18  PM 

TIN  File 

APP-TIN-2001- 

CORR.NPT 

07/03/14 
10:53  AM 

APP-TIN-2003- 

CORR.NPT 

07/03/14 
10:52  AM 

APP-TIN-2010- 

CORR.NPT 

07/03/14 
10:54  AM 

APP-BR2- 

TIN.NPT 

02/24/14 
3:35  PM 

APP-BR2- 

TIN.NPT 

02/24/14 
3:35  PM 

APP-BR2- 

TIN.NPT 

02/24/14 
3:35  PM 

APP-BR3- 

TIN.NPT 

02/24/14 
3:35  PM 

APP-BR3- 

TIN.NPT 

02/24/14 
3:35  PM 

APP-BR3- 

TIN.NPT 

02/24/14 
3:35  PM 

QOT  File 

APP-QOUT- 
2001. NPT 

02/25/14 
1:33  PM 

APP-QOUT- 
2003. NPT 

02/25/14 
1:33  PM 

APP-QOUT- 
2010. NPT 

05/25/14 
1:33  PM 

QDT  File 

APP-QDT-2 001- 
ADJUST.  NPT 

12/23/14 
2:11  PM 

APP-QDT-2003- 

ADJUST.NPT 

01/14/15 
2:10  PM 

APP-QDT-2  010- 
ADJUST.  NPT 

01/20/15 
10:29  PM 

TDT  File 

APP-TDT- 
2001. NPT 

07/03/14 
10:53  AM 

APP-TDT- 
2003. NPT 

07/03/14 
10:52  AM 

APP-TDT- 
2010. NPT 

07/03/14 
10:54  AM 
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Table  B3.  Inventory  of  files  needed  to  run  the  APPLPM  (Predictive  Model). 


Run  Name 

CYOl-USGS-PortRunlO 

CY03-USGS-PortRunl3 

CYlO-USGS-PortRu  n07 

File  Type 

Calibration  - 
2001 

Date 

Stamp 

Verification  - 
2003 

Date 

Stamp 

Verification  - 
2010 

Date 

Stamp 

W2_C0N.NPT 

“ 

10/26/15 
2:10  PM 

- 

01/09/15 
10:57  AM 

“ 

10/26/15 
3:57  pm 

QOT  File 

APP-QOUT- 

2001- 

PortRunOl.NPT 

05/05/15 
9:41  AM 

APP-QOUT- 

2003- 

PortRunOl.NPT 

APP-QOUT- 

2010- 

PortRunOl.NPT 

W2_SELECTIVE.NPT 

-- 

08/03/15 
2:10  PM 

-- 

08/03/15 
2:10  PM 

-- 

08/03/15 
2:10  PM 

DYNSPLIT_SELECTIVE1.NPT 

-- 

05/05/15 
1:37  PM 

-- 

05/05/15 
1:37  PM 

-- 

05/05/15 
1:37  PM 

DYNSPLIT_SELECTIVE2.NPT 

- 

05/05/15 
1:37  PM 

- 

05/05/15 
1:37  PM 

-- 

05/05/15 
1:37  PM 

DYNSPLIT_SELECTIVE3.NPT 

- 

05/05/15 
1:37  PM 

- 

05/05/15 
1:38  PM 

-- 

05/05/15 
1:37  PM 

**Note:  The  same  w2_selective.npt  AND  dynsplit*  files  are  used  for  all  3  cases.  Unless  noted  above,  the  files  used  in  Table  B2 
apply  to  the  Predictive  Model. 
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Appendix  C:  APPLM  and  APPLPM  Files 

This  appendix  serves  to  provide  a  description  of  each  file  needed  to  run  the 
model.  The  files  are  grouped  by  year.  The  ERDC  typically  has  the  following 
file  organization  system  (see  Table  Cl). 


Table  Cl.  Typical  file  organization. 


CY01 

Main  folder  for  year  identification  for  the  particular  model.  Most 
models  will  be  designed  to  run  with  multiple  years. 

Results 

Upon  running  the  model,  the  results  are  moved  out  of  the 
executables  folder  and  into  their  own  folder;  typically,  these 
folders  are  named  something  like  CYXX_RunXX.  NOTE:  Always 
copy  the  control  file  used  for  the  run  into  the  results  folder  so  that 
you  can  duplicate  the  run  in  the  future  if  necessary. 

Executables 

This  is  where  all  of  the  necessary  files  needed  to  run  the  model 
are  located:  W2  executables,  Inflows,  Outflows, 
Temperature/Concentration  files,  Met  files,  Bathymetry,  etc. 

Table  C2.  Files  needed  to  run  APPL  Model  for  each  year. 


File  Description 

CY01 

CY03 

CY10 

Graph  File 

graph,  npt 

graph,  npt 

graph. npt 

Control  File 

w2_con.npt 

w2_con.npt 

w2_con.npt 

Bathymetry  File 

APP-BATH-FINAL- 
2001. NPT 

APP-BATH-FINAL- 
2003. NPT 

APP-BATH-FINAL- 
2010. NPT 

Meteorology  File 

APP-M  ET-2001.  N  PT 

APP-M  ET-2003.  NPT 

APP-M  ET-2010.  NPT 

Wind  Sheltering  Coefficient  File 

APP-WSC.NPT 

APP-WSC.NPT 

APP-WSC.NPT 

Shade  File 

APP-SHD-2.NPT 

APP-SHD-2.NPT 

APP-SHD-2.NPT 

Upstream  Inflow  File 

APP-QIN-2001.NPT 

APP-QIN-fromGates- 
2003. NPT 

APP-QIN-fromGates- 
2010. NPT 

Upstream  Temperature  File 

APP-TIN-2001- 

CORR.NPT 

APP-TIN-2003-C0RR.NPT 

APP-TIN-2010-C0RR.NPT 

Branch  2  Inflow  File  (zero) 

APP-BR2-QIN.NPT 

APP-BR2-QIN.NPT 

APP-BR2-QIN.NPT 

Branch  2  Temperature  File 
(placeholder) 

APP-BR2-TIN.NPT 

APP-BR2-TIN.NPT 

APP-BR2-TIN.NPT 

Dam  Outflow  File 

APP-Q0UT-2001.NPT 

APP-Q0UT-2003.NPT 

APP-Q0UT-2010.NPT 

Distributed  Tributary  Inflow  File 

APP-QDT-2001- 

ADJUST.NPT 

APP-QDT-2003- 

ADJUST.NPT 

AP  P-QDT-2  010- 
ADJUST.  NPT 

Distributed  Tributary 

Temperature  File  (duplicated 
upstream  temps) 

APP-TDT-2001.NPT 

APP-TDT-2003.NPT 

APP-TDT-2010.NPT 
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